Poincare recurrences and transient chaos in systems with leaks 
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In order to simulate observational and experimental situations, we consider a leak in the phase 
space of a chaotic dynamical system. We obtain an expression for the escape rate of the survival 
probability applying the theory of transient chaos. This expression improves previous estimates 
based on the properties of the closed system and explains dependencies on the position and size of 
the leak and on the initial ensemble. With a subtle choice of the initial ensemble, we obtain an 
equivalence to the classical problem of Poincare recurrences in closed systems, which is treated in 
the same framework. Finally, we show how our results apply to weakly chaotic systems and justify 
a split of the invariant saddle in hyperbolic and nonhyperbolic components, related, respectively, to 
the intermediate exponential and asymptotic power-law decays of the survival probability. 



I. INTRODUCTION 



The numerical and theoretical study of dynamical sys- 
tems can resemble observational and experimental situ- 
ations if the escape of an initial ensemble of trajecto- 
ries through a leak placed in the otherwise closed phase 
space is considered. The idea of introducing a leak in 
a closed chaotic dynamical system to generate transient 
chaos was first suggested by Pianigiani and Yorke 
Later this problem was discussed in detail in the context 
of fractal exit boundaries @ , of geometrical acoustics [|[ , 
of quantum chaos of controlling chaos [j|, of reset- 
ting in hydrodynamical flows M, 0] , of leaked Hamilto- 
nian systems Q , of astronomy [9J and of cosmology [Icl ] . 
The subject has recently received a renewed attention 
[O, El 111 EI EE EE E3- Part of this renewed inter- 
est comes from recent developments of quantum chaos, 
that started to consider carefully the effect of measure- 
ment devices, absorption, and other forms of leaking both 
theoretically and experimentally [ill El El, EU • It is by 
now clear that leaking dynamical systems provides a way 
to 'peeping at chaos' Jla|. and can be used as a kind of 
chaotic spectroscopy j ij . 

From a theoretical point of view, leaking systems pro- 
vide an interesting bridge between open and closed dy- 
namical systems. In closed systems the phase space co- 
ordinates remain confined, and a rigorous mathematical 
approach based on attractors and ergodic components 
can be employed [Lsj . An important quantity is the dis- 
tribution _£r (T) of the first Poincare recurrence times T 
HE EEL Hi, Hi HH to a preselected region of the 
phase space. As originally proposed by Chirikov and Shc- 
pelyansky [25[ (see also [23|), p r (T) is a useful analyser of 
the entire dynamics. In contrast, in open systems trajec- 
tories may leave the region of interest (e.g., tend to ±oo) 
and the method of transient chaos is usually employed 
[l8l |26| . In this case, the dynamics is characterized by 
the escape time distribution p e (n) (the derivative of the 
survival probability) . 

Chaotic dynamics usually leads to an asymptotic ex- 



ponential decay of both p r and p e 



Pr,e(f) 



-tr.et 



(1) 



for large enough t — {T,n}. Exponentials are the sig- 
nature of strong chaotic properties as seen, e.g., in the 
divergence of nearby initial conditions, in the decay of 
correlations, and in the convergence to equilibrium dis- 
tributions. The exponent 7^ is the key quantifier of 
system specific characteristics and will be considered in 
detail in this paper. 

The possibility of comparing a system with leak to its 
corresponding closed system shows numerous advantages 
in comparison to naturally opened systems. For instance, 
the rate 7 rje in Eq. |T]) can be estimated from the prob- 
ability of escaping or returning. For small leaks, this 
probability is obtained by computing the closed systems 
natural measure \i of a typical leak region /, and it follows 

that [Emm 



7 r = 7 e = fl(I) 



1 



for 1 > fj,(I) > 0, (2) 



where (t) denotes the mean recurrence or escape times, 
(T) or (n). 

Despite its simplicity, relation @ is of little practical 
relevance since the limit fj,(I) — > is never achieved in ei- 
ther numerical or experimental applications [H IH, Hi| . 
Here we do not restrict ourselves to this limit and find 
that, apart from being unrealistic, it masks different in- 
teresting relations. We perform a formal description 
through the theory of transient chaos, i.e., in terms of 
a nonattracting chaotic set [l||] and the conditionally in- 
variant measure P,H3|, and obtain distinct relations for 7 
and 1 / (t) as a function of the position and size of the re- 
gion / and of the initial ensemble. Relation © does 
not hold, but we show that 7 r = 7 e = 7 remains valid. 
This is based on a full correspondence with the recur- 
rence problem obtained for a special initial ensemble and 
setting the leak to be equal to the recurrence region of the 
closed system. We have explored the most striking conse- 
quences of this unified treatment from the point of view 
of Poincare recurrences in a previous short paper [3l| . 



2 



Here, after describing the relationship between the two 
problems with additional details, we focus on the escape 
problem and consider the recurrence problem as a par- 
ticular, yet important, case. We also use this connection 
to apply and adapt a method to determine 7 [27J based 
only on the escape probabilities for short times. This 
is a main advantage over genuine open systems where 
long-term periodic orbits are needed to determine j e . 

In systems showing weak chaos, as e.g., Hamiltonian 
systems with mixed phase space, the exponential decay 
law ([T]) experiences a cross-over for longer times towards 
an asymptotic power-law behavior [23l. [24. l32| . We argue, 
however, that on intermediate times a decay rate 7 can 
be well defined 



Pr, e {t) (X 



e 



7* 



for intermediate t < t c 
for t > t r . 



(3) 



We interpret our results based on an effective splitting 
of the chaotic saddle into hyperbolic and nonhyperbolic 
components. 

The paper is organized as follows. In Sec.|TT]we review 
the main properties of Poincare recurrence times. The 
corresponding results for escape times appear in Sec. IIIIl 
together with the dependence of 7 and (t) on /. The 
equivalence between recurrence and escape problems is 
carefully discussed in Sec. |IVl In Sec. [V] numerical il- 
lustrations of the results are provided for paradigmatic 
dissipative models. In Sec. I VII the nonhyperbolic case is 
investigated using an 

area-preserving map. Relations to periodic orbits and 
an efficient method to determine 7 are described in 
Sec. IVII1 Finally, conclusions and discussions appear in 
Sec. IVIII1 where the case of higher order recurrences is 
also investigated. 



II. POINCARE RECURRENCES IN CLOSED 
SYSTEMS 

Poincare recurrences to a certain region of the phase 
space have played an important role since the founda- 
tion of the kinetic description of noneqilibrium processes 
[2lL |24| ]. We define the problem in the context of low- 
dimensional chaotic maps that can be dissipative [33j or 
Hamiltonian [13, Hi| . 

Consider a discrete time chaotic system x n +i = M(x n ) 
defined on a bounded phase space x G T. Let the natural 
ergodic measure on the invariant chaotic set (e.g., strange 
attractor or chaotic sea) be denoted by p and its density 
by pp. We define the recurrence region as a subset I C 
r with p(I) > 0. As stated above, we do not restrict 
ourselves to the unrealistic limit p(I) — > [27L I2H . [2^ |. 
The Poincare recurrence theorem ensures that for almost 
all initial conditions xq G /, there are infinitely many 
time instants n = n\,n2, ■ ■ ■ such that M n (xo) G The 
first recurrence times are defined as Tj = rii — nj_i, with 
i > 1 and no = (if the points remains in / at the ith 
iterate Tj = 1). The recurrence time distribution p r {T), 



T > 1 , is the probability of finding Tj = T in an infinitely 
long trajectory. The cumulative version can be written 
as Pr(i") = Y^t=t P r (-^0 • The mean recurrence time is 
defined as 



(T) 



lim 

N^oo N 



1 JV 00 



1=1 



T=l 



Due to the ergodicity of the measure p, the specific choice 
of the initial point xq G I is irrelevant for p r (T). In- 
stead of following a single trajectory, p r {T) can also be 
obtained by using an ensemble of trajectories chosen in- 
side / according to the invariant density p^. A conse- 
quence of ergodicity is Kac's lemma 19]: 



(T) 



1 



(4) 



which is valid for any recurrence region I. 

The recurrence time distribution of an uncorrelated 
random process with fixed return probability p follows 
a binomial distribution (27j 

p r (T)=Ml-/i) T - 1 = T 4-^ (1 ^ )T (5) 
1 p 

The relaxation rate 7 r defined in Eq. {I} is thus given for 
such a random process by 



7,; = -ln(l - p). 



(6) 



For small recurrence probability p — * 0, 7* — > p, and the 
Poisson distribution p r {T) — penp(-pT) is obtained. 
Identifying p with p(I) (recurrence probability equals to 
the natural measure of the recurrence region), we see 
that a Poisson process is the simplest process leading to 
Eq. @. Indeed, the Poisson distribution is proved to 
describe deterministic hyperbolic systems (see references 
in gll). 

We study next the distribution of recurrence times 
p r {T) for generic chaotic systems. In agreement with 
Eq. |T]), we write 



Pr(T) 



irregular 

9re~ lrT 



for 1< T < T* 
for T > T*. 



(7) 



Implicit is the view that the approach toward the expo- 
nential distribution for T — > 00 does not occur uniformly, 
but that after some short time T* the distribution is 
practically very close to an exponential [27]. This will 
be theoretically justified in the next section and numeri- 
cally illustrated in Sec.[V] We are interested in the decay 
rate j r , which typically deviates from the random esti- 
mate ([6]) , and also in the short time irregular fluctuations 
for T < T*, where T* is of the order of (T) and will be 
better interpreted in Sec. |lVl These short time irregular 
oscillations [denoted as irregular in Eq. ([7]) and equations 
hereafter] depend strongly on the choice of the recurrence 
region /. For instance, if a periodic orbit of (short) pe- 
riod p is present inside /, there is a higher probability of 



3 



recurrence at T = p. As shown in Ref. [27] and will be 
discussed in Sec. I VIII normalization and Eq. ^ imply 
that the short-time fluctuations uniquely determine the 
deviation of j r from 7*. Therefore, for different recur- 
rence regions I of the same natural measure p(I) [while 
(T) is fixed and given by (fflj, the value of j r depends 
strongly on the position of I (27| . 



III. ESCAPE IN OPEN SYSTEMS 
A. General relations for transiently chaotic systems 

By now, transient chaos has a well established theory 
In this case, the discrete system x n +i — M(x n ) 
is defined in an unbounded phase space T. A region of 
interest icT (non-trivial dynamics) is defined in such 
a way that trajectories may leave A \M(A) ^ A] but do 
not return to it [M(T \ A) F\A = 0]. An initial ensemble 
is started according to a smooth density po(x), with x £ 
A. The escape time distribution p e (ri) is given by the 
fraction of trajectories that leave A at time n > 1. For 
chaotic systems it can be written as 



, x _ J irregular for 1 < 
PeW - I g ee -i.n for n > 



1 < n < n*, 
n> n* 



(8) 



where n* corresponds to a short convergence time and 7 e 
is the escape rate. The survival probability inside A up 
to time n is given by P e (n) — Y^n+i Pe{ n ')-> an( A is also an 
exponential distribution with the same j e as in Eq. l[5]). 
The average escape time, (n) e = n Pe( n ), also called 
the lifetime of chaos, is usually estimated by the recipro- 
cal of 7 e [H, |2f|. Note, however, that the exact average 
is, in general, different from l/7 e and depends (just like 
n* e ) on p . 

The ergodic theory of transient chaos explains the 
above properties by the existence of a non-attracting 
chaotic set. In invertible systems this set is, e.g., a chaotic 
saddle [56[: the set of orbits that never escape either for- 
ward or backward in time 

[HIHHI]. For systems with 
strong chaos, the saddle is an invariant fractal set of zero 
measure. Long-time transients correspond to trajecto- 
ries that approach it closely, along its stable manifold, 
and leave it along its unstable manifold. Being an invari- 
ant object, the chaotic saddle is, of course, independent 
of the initial density pg. Therefore, assuming that the 
distribution po is non-vanishing in at least some region 
around the stable manifold of the saddle (and that there 
are no disjoint saddles) the escape rate 7 e is independent 
of the initial density. It can be related to the geometrical 
and dynamical properties of the saddle by the Kantz- 
Grassberger formulas [3~ij , which for two-dimensional in- 
vertible maps are 



7 e = A(l - d u ), Xd u = ~X'd s , 



(9) 



where A (A') is the positive (negative) Lyapunov expo- 
nent on the saddle and d u (d s ) is the information dimen- 
sion of the saddle along its unstable (stable) manifold 



[la . |26[ . The time n* , below which orbits escape without 
approaching the saddle, can be estimated by the inverse 
of the negative Lyapunov exponent A'. 

A rigorous description of open systems is possible in 
terms of the conditionally invariant measure P, l30l . l35l . 
HU, hereafter called c-measure. A measure p c is said to 
be conditionally invariant if 



Pc(M~\E)) 
p c (M-i(A)) 



Pc(E), 



(10) 



for any E C A. In words: the c-measure is not invari- 
ant under the map: p c (Al^ 1 (E)) 7^ p c (E), but is pre- 
served under the incorporation of the compensation fac- 
tor in the denominator of (jTUJ) . The compensation factor 
p c (M~ 1 (A)) < 1 corresponds to the c-measure of the 
set remaining inside A, i.e., the trajectories that do not 
escape in one iteration. The c-measure of the region of 
interest is p c {A) = 1. The c-measure has a density p c , 
which is unique and the only attractor for smooth initial 
densities p . Taking into account the exponential escape 
expressed in (JHJ, one sees that convergence is only possi- 
ble if Q 



le = -\n(p c (M-\A))). 



(11) 



Considering the evolution of po through an evolution op- 
erator, the Perron-Frobenius operator 0j[13]> ex P(7e) is 
the largest eigenvalue and p c is the corresponding (right) 
eigenvector. The c-measure describes the escape process 
of the system and therefore p c is concentrated on the 
saddle and along its unstable manifold [3(3]. 

Numerically, p c is obtained multiplying pa at each it- 
eration by a constant factor which, after many itera- 
tion, is necessarily equal to exp(7 c ). In practice, we it- 
erate a large number of trajectories and we obtain p c by 
renormalizing the surviving trajectories at a long time. 
p c (M~ 1 (A)) in (jnj [or 1 - p c (I) in HU) below] is ob- 
tained as the fraction of all surviving trajectories that 
escape in the next time step or, equivalently, by fitting 
7e to p e (n) and inverting (fTTj) . 



B. Systems with leaks 

Open systems M can be obtained from the closed sys- 
tem M defined on T, described in Sec. [TH by introducing 
a leak in a region I CT 



M(x n ) 



M(x n ) 
escape 



if x n £ I 
if x n £ L 



(12) 



Notice that, since the escape happens one step after en- 
tering /, the map M is defined in / and initial conditions 
can be in /. This procedure defines an open system as 
that of Sec. HTlAl with A = T. 

An estimate 7* of the escape rate 7 e in systems with 
leaks is usually given in terms of the Frobenius-Perron 
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relation [f| [2(| which expresses that the number of par- 
ticles not escaping within one time step is proportional 
to the natural measure outside of the leak: exp (—7*) = 
1 - fi(I), i.e., 



7e 



ln(l- M (/)) [«M0Ofor MW-0], (13) 



where p(I) is the natural measure of leak /. This re- 
lation is equivalent to the binomial estimate for recur- 
rence times ([H]) and we call it hereafter as the naive es- 
timate of 7 e . Deviations of 7 e from 7* were reported in 
Refs. [U, 0] and associated to the existence of short-time 
periodic orbits of the closed system inside /. 

The results of Sec. lIII Al on open systems can be applied 
to systems with leaks. Since M^ 1 (A) = A\ I, it follows 
that ^ C (M _1 (A)) = fx c (A) - p, c (I) = 1 - fi c (I), and the 
following simple relation is obtained from pip 



7e = -ln(l-M c (7)). 



(14) 



This exact expression, which was previously obtained in 
Ref. 0], corresponds to the naive estimate (fT3"|) replac- 
ing p(I) by /J-c(I)- They become equivalent in the limit 
of small leak region p(I) — * because p c tends to p 
[37L [381 . It is the c- measure of the leak which determines 
the escape rate. Even though numerically such a mea- 
sure is easily calculated. Proofing the existence of such 
conditionally invariant measure may be a very involving 
mathematical task 0, [35| (see, e.g., Ref. [36| for a class 
of maps with leaks). Here we employ a pragmatic ap- 
proach and implicitly assume that such a measure exists 
in the large class of systems with leaks where numerical 
simulations indicate exponential decay of p e (n). 



C. Dependence of the mean escape time on the 
initial density 

While the escape rate 7 e is independent of the initial 
density po, as expressed by Eq. (TTJ}, the mean escape 
time (n) e and the full distribution p e (n) change consid- 
erably with pq. It is thus instructive to discuss the fol- 
lowing special initial conditions: 

Conditionally invariant density po = p c : from a for- 
mal mathematical perspective the natural choice of ini- 
tial conditions is according to the density of the c- 
measure p c (x), which can be considered as the natural 
measure of the open system. This corresponds to setting 
the initial density in agreement with the escape process. 
Therefore, no oscillations are observed in Eq. ([8]), i.e., 
n* = 1. Since 7 e is given by (fT4"]) and the distribution is 
normalizable, p e (n) follows a binomial distribution 



Pe(n) = A*c(l - Mc)" : ' 



1 -Mc 



with p c = p c {I), in analogy with Eq. ([5]). The mean 
escape time is in this case 



{n) c 



1 



1 



p c {I) 1 - e-i- 



with po = Pc- (15) 



Natural density po — p^: coherent with the idea of 
systems with leaks, the initial density is chosen accord- 
ing to the natural measure of the closed system over the 
full r. It corresponds to leaking the system after equilib- 
rium has been achieved in the closed system. We show 
in Sec.|V] that 

( n )v ~ TTTn = 1 with ' ^0 = /V ( 16 ) 



A*c(J) 



1 - e-T« 



Similar results are obtained for densities proportional 
to p^ (equal to p M apart from normalization in a sub- 
space of r). 

Smooth density po = p s : in some systems, and for 
numerical simulations, it is natural to define the initial 
density to be a smooth function over T or a part of T. 
A particular case is that of a uniform density. Estimate 
(TTB]) remains valid for such initial densities as well, as 
shown in Sec. IVl 



IV. RELATION BETWEEN RECURRENCE 
AND ESCAPE 

In this section we show that there is an initial density 
Po = Pr for which p r {T) = p e {n) for T = n. We use this 
special initial condition to establish a relation between 
recurrence and escape. 



M(I)= escape 




FIG. 1: (Color online) Illustration of the equivalence between 
recurrence and escape times. 

Density equivalent to recurrence po — p r : consider that 
trajectories are injected from outside through the leak. 
More precisely, consider the distribution p r obtained as 
the first iterate through map M{x) of the points x E I 
distributed according to the natural density p^. This can 
be obtained applying the Perron-Frobenius operator [13] 
as 

where M~ 1 {x)C\I denotes the points that come from /, J 
is the Jacobian of the map, and p(I) ensures normaliza- 
tion (5?| . In the case of invertible maps, A/ -1 is unique 
and the density (JTTj) is equivalent to 



p r (x) 




if xeM(T), 
else. 



(18) 
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Large leaks 


Small leaks 


TVT P fl c; 11 r P 

IVIG CXO LI! t/O 


Finite n(T} 11 ( T\ 


u(t\ — n (r) — > n 


Escape rate 


T = -In (1 - /!„(/)) / - ln(l - /!(/)) = 7 * 


7 = p c (/) = p(J) 


Mean time 
Po 


{n)r = p(I) * 1 - e-7 (n)c = ^(J) = 1 - e- (n) - s * p c (J) = 1 - 
Recurrence: p r c-measure: p c natural, smooth: p^s 


/ \ 1 1 

( " > = K0 = 7 



TABLE I: (Color Online) The escape/recurrence rate 7 = 7 r = 7 e and the mean escape time (recurrence time for initial 
density p r ) (n) expressed in terms of the natural invariant measure p(I) and the conditionally invariant measure p c (I) of the 
leak/recurrence region /, for large and small leaks. 



Consider now an infinitely long trajectory used to calcu- 
late the recurrence times T described in Sec. |TTJ In view 
of the ergodic theorem (time average equals ensemble av- 
erage) , we see that one iteration after returning to / the 
points of this trajectory are distributed precisely as p r . 
As illustrated in Fig. [TJ due to the one step time shift 
inserted in definition (|12p , all escape times n of p r corre- 
spond to a recurrence time T[58|. Since this is valid for 
all n, p e (n) = p r (T) for p = p r . In particular, j r = j e , 
and since 7 e is independent of the initial density, this im- 
plies that for a fixed recurrence/escape region / the decay 
rate of the Poincare recurrences and the escape rate of 
the corresponding leaked system coincide: 



7r = 7e 



7- 



(19) 



The mean recurrence time is given, however, according 
to Kac's lemma (@| as 



(T) 



with po = p r . 



(20) 



The mean recurrence time is thus determined by the nat- 
ural measure p(I) and deviates considerably from the 
typical mean escape time given by (1161) . The reason is 
that p r in 1|17|) is very atypical from the point of view of 
the c-measure. In fact, p c is concentrated along the sad- 
dle's unstable manifold, i.e., points x such that M~' l (x) 
never escape. On the other hand, all points on p r came 
from the leak and therefore the support of p r does not 
overlap with the unstable manifold of the saddle (but it 
does with the stable manifold). 

In view of the results above, we consider hereafter 
the recurrence problem as an escape problem with p = 
p r [59l | . Coherently, we omit the indices r and e, and T 
will also be denoted by n. The values of the escape rate 7 
and of the mean escape time (n) for the different initial 
densities p are summarized in Table HI 

It is now possible to compare the results of Ref. [5( for 
escapes with those of Ref. [27j for recurrence. In both 
cases the fully chaotic logistic map was carefully investi- 
gated and deviations of 7 from the naive estimate 7* were 
reported. The deviations were associated with the pres- 
ence of periodic orbits inside I and interpreted in terms 
of the overlap of the pre- images of / (in Ref. ||) and 
of the short time oscillations (in Ref. [23] )• In our uni- 
fied perspective, results of Refs. 27 1 and Q correspond 



to choosing initial densities p r and p^ or p s , respectively 
(see Table |T|. Since the escape rate is independent of this 
choice, we see that both references are giving intuitive in- 
terpretations for the deviation of the actual 7 from the 
naive estimate 7*. Analogous observations in Hamilto- 
nian systems have been reported in (recurrence) and 
in [tJ (escape). In a different approach, 7 has been re- 
lated to the diffusion coefficient in both recurrence (23j 
and escape [39j problems. 

Formally the problem can be described in terms of the 
c-measure and the convergence to it. The difference be- 
tween 7 and 7* correspond to the difference between p(I) 
and p c (I), as given by Eqs. (flU)) and (fT¥jl . As empha- 
sized in our previous publication [31] , this is surprising 
since a measure of open systems is employed to describe 
properties of closed systems (recurrence). In practice, 
initial densities are most often taken proportionally to 
the natural density p^ or to a smooth density p s . In such 
cases we expect a fast convergence to p c , and this is why 
l/Hc{I) [Eq. (US}] is a better estimate of (n) [in Eq. (fLo|) ] 
than l/p(I) [Eq. ([3}, which is obtained by the atypi- 
cal po = p r ]- The escape rate is related to p c (I) which 
is itself specially sensitive to the choice of the location, 
shape, and size of the finite leak I. It is natural to expect 
that p c {I) < p{I) since the saddle exists only in A\I and 
therefore p c is concentrated only along its unstable man- 
ifold. This implies that the mean recurrence time (T), 
i.e., (n) with po — Pr should be shorter than the mean 
escape time (n) for a typical po- This statement is rigor- 
ously proved in Ref. [40j for po = Pn- 



V. NUMERICAL RESULTS 

We illustrate the previous results through numerical 
simulations of the dissipative Henon map 



x n+ i = 1 - \Ax Ti 



y n+ i = 0.3x n . 



(21) 



Initial conditions inside a large basin of attraction close to 
the origin converge to the well studied Henon attractor, 
shown in Fig. [DJa). We let trajectories leak the system 
through the leak region /: circle of radius S centered at 
some (x c ,y c ) on the attractor. The Sinai-Ruelle-Bowen 
measure is the natural measure p describing the closed 
system. All the initial densities discussed in Sec. IIII CI 



FIG. 2: (Color online) Henon map (I21|l with a leak 
at (x c ,y c ) = (1.0804,0.0916) of radius 5 = 0.1, depicted 
as a red circle, (a) Different initial densities: p^ of the 
natural measure on the Henon attractor, p s — const, on 
< x < 1, —0.2 < y < 0.2, and p r , equivalent to recurrence. 
The inset shows the nontrivial dependence of the projection 
p r (y) of p r on the y axis, (b) Invariant saddle composed of 
the points that do not enter I for forward and backward iter- 
ations, (c) Unstable manifold of the invariant saddle (p c ) and 
(d) stable manifold in grey (orange) of the invariant saddle. 



are marked in Fig. [5J Panels (b)-(d) show the invariant 
saddle and its manifolds. A new feature in comparison 
with saddles of naturally open systems is that the unsta- 
ble manifold is not a single fractal curve but consists of 
disjoint pieces cut by the leak and its images. 

In Fig. [3] we show the escape time distribution for the 
initial densities depicted in Fig. [21 It is apparent that all 
initial densities lead to different short time behavior but 
to the same escape rate 7, which is significantly different 
from the naive estimate (fT3l). 



We investigate now the dependence of 7 and p(I) on 
the radius 5 of the leak. Since the natural and the c- 
measure have different fractal properties, for small (but 
not infinitesimally small) 5 one expects the scaling rela- 
tion 



1 



5 D \ 



(22) 



where D\ is the information dimension of the chaotic at- 
tractor, d s is the partial information dimension of the 
saddle (system with leak) along the stable direction (the 
c measure is smooth along the unstable one), and © 
has been used. The numerical comparison can be seen 
in Fig. [4] where the inverse mean escape time is plot- 
ted for different initial densities. For the c-measure and 
S <C 1, corresponds to p, c (I) through relation ([T4"|l . 

while l/(n) = p(I) with p — p r . The deviation between 



FIG. 3: (Color online) (a) Escape time distributions for the 
Henon map (|21[) with a leak of radius S — 0.05 at (x c ,y c ) = 
(1.0804,0.0916). Results for initial densities Pr,p M > and p s 
(shown in Fig. [2} are presented from top to bottom, respec- 
tively. The thin solid line corresponds to 7* = 0.0385 given 
by ((6J and strongly differs from the actual exponent 7 = 
0.545. The inset shows the short time behavior, (b) The 
same as in (a) multiplied by exp(7n). Three different pro- 
cedures were employed to fit 7 using p(n) for n < n* (see 
Sec. IVIl)) and are indicated by the dashed lines (i)-(iii). In all 
cases n* = 35 and (i) corresponds to the use of relation (|33|l . 
(ii) follows from relations (|31[l and f|32p with the numerical 
value of (n) for po — p r , and (iii) from relations (|31[) and (|32[) 
with the numerical value of (n) for po = Ps- 



these two curves is a measure of the error of the naive 
estimate (fl~3j) . The general dependence, for 5 — ^ 0, is 
in agreement with the expected 5~ Dl relation, by taking 
into account that Di ;$ Do, where Do is the fractal di- 
mension of the attractor. From the inset it is apparent 
that initial conditions and p s according to the natural 
measure and to the homogeneous distribution, respec- 
tively are closer to the results obtained with the density 
p c of the c-measure than that obtained with initial den- 
sity p r of the recurrence problem. In other words, the 
mean escape time for typical smooth densities is better 
approximated by Eq. (JT6J) than by l/p(T) [Eq. (gj)]. 

A systematic investigation of the dependence of 7 on 
the position of the leak in the Henon map would be very 
involved. Apart from the difficulty of having two dimen- 
sions of the phase space, the natural (SRB) measure of 
the system is not known a priori. Therefore, in order 
to fix the natural measure of the leak p(I), a different 
value of S should be chosen for each new (x c , y c ). We in- 
vestigate therefore this issue in the fully chaotic logistic 
map 



Xn-\-i — 4x n (l 



(23) 



In this case the density p M of the natural measure of the 
system is given by [18[ 



P ,{X)='K- 1 [X{1-X)\- 1 ' 2 . 



(24) 
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FIG. 4: (Color online) Dependence of the inverse mean es- 
cape time l/(ra) on the radius S of the leak for fixed (x c , y c ) — 
(0.015714,0.268892). Initial densities equivalent to those 
shown in Fig. [2] are taken. The straight line corresponds to 
<5~°°, where Do = 1.27 is the fractal dimension of the attrac- 
tor [53j] . The inset shows a magnification of the main graph 
for large 8. 

The results are shown in Fig. [Hb) and confirm the strong 
dependence of 7 on /, as well as the conclusion drawn 
from Fig. [4] that for general smooth initial densities 7 is 
better approximated by (|16p than by l//x(J), reflecting 
once more the fast convergence of po to p c . Moreover, 
it is clear that for most x c one has (n) e > (n) r , i.e., 
/i(J) > p c {I), in agreement with the discussion at the end 
of Sec. IIVI The remarkable exceptions are related to the 
existence of short time periodic orbits inside the leak, as 
previously described in Refs. d, Ezl Ezl • An explanation 
for these cases can be given as follows (see Sec. IVIII for 
details): if a periodic orbit of period n p is inside / we 
expect a high recurrence probability P(T = n p ); since 
(T) is fixed by (U|), we predict that 7 < 7* and therefore 
(n r )/(n e ) < 1, as observed in Fig. OJb). 




- , . 1 . p Q =p r : recurrence, 1/<n> = |i(l) 

|| 0.8" 2. p =p : c-measure, 1/<n>= u (l)=1-e^ 

0) 

A " □ □ 3. p =p : inv. measure 



V 0.7 * * 4- p =p s : uniform in [0,1] 

a" 
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x 
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FIG. 5: (Color online) (a) Graphical representation of the 
logistic map (123 p with a leak I at x c — 0.29427 of half-width 
(5 = 0.00716 so that p(I) = 0.01. The dashed-line corresponds 
to Pfj,(x), Eq. (I24|l . and the (blue) histogram to a numerical 
distribution of the c-measure of the chaotic repeller obtained 
by leaking the system at I. The dotted (red) histogram cor- 
responds to p r (x) given by (|17p multiplied by p(M(I)). (b) 
Dependence of the inverse mean escape time, multiplied by 
{n) r , on the leak position x c in the logistic map ([23} at an in- 
terval length such that p(I) — 0.01. The initial densities used 
are (see Table UJ: p r : equivalent to recurrence, p M : natural 
density Eq. (]24p . and p s : chosen uniform in < x < 1. 



VI. NONHYPERBOLIC HAMILTONIAN 
SYSTEMS WITH POWER-LAW TAILS 

A. Theory 

Our discussion has covered so far the case of strongly 
chaotic systems where asymptotically in time, the dis- 
tribution of escape and recurrence times is exponential 
(for certain classes of systems this has been formally 
proved (44|). We concentrate in this section on weakly 
chaotic systems where deviations from exponential de- 
cay are observed for long times. The simplest exam- 
ples of weakly chaotic systems are one-dimensional in- 
termittent maps with a marginal fixed point, where an 
asymptotic power-law decay is proved [45l |46| . The most 
prominent examples are, however, Hamiltonian systems 
with mixed phase space, where asymptotically a power- 
law decay is observed. This occurs in systems possessing 
isolated marginally unstable orbits [471 ] . sharply divided 
phase space [481 ] . hierarchical Kolmogorov-Arnold-Moser 



(KAM) islands coexisting in a chaotic sea pll. I2H H2I [43j . 
and also in higher dimensions 49] . The power-law behav- 
ior is due to the nonhyperbolicity of the dynamics and it 
is present both for recurrences and escapes. The power- 
law exponent is independent of the choice of / and po, 
provided the latter is concentrated away from the non- 
hyperbolic regions. Below we show that the same holds 
for the intermediate time exponential decay of the recur- 
rence and escape time distributions, for a given I. 

The distribution of escape and recurrence times in 
Hamiltonian systems with mixed phase space can thus 
be given as 

{irregular for n < n*, 

jae~ jn for n* < n < n a , (25) 

7[ae~ 7 ™ + b(~fn)~ a ] for n > n a , 

where ae _7 ™ a 3> b(jn)~ na . The factor 7 is written out 
on the right hand side in order to facilitate the connec- 
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tion to the continuous time limit. Contrary to what has 
previously been claimed [50], 7(0 + b) ^ 1 in Eq. ([25]) . 

The question whether the asymptotic regime has a well 
defined (universal) power-law is still under investigation 
for the case of area-preserving maps (see Ref. [43j for 
the latest result that indicates that a ~ 2.57). Here, for 
simplicity we write the asymptotic decay as n~ a , but it 
is meant to describe the power-law like behavior usually 
observed. 

We define the cross-over time n c between exponential 
and power-law decay as: 

ae-^ = b(~fn c )- a p(n c ) = (26) 

In the framework of recurrences, the different times in 
Eqs. (|25|) and l[26|) can be interpreted as: 

n* : is the short-time memory, within which fluctuations 
due to short time periodic orbits appear. After this 
time the recurrences loose correlation and a ran- 
dom approximation for the return is valid (Poisson 
process). 

n a : is the minimum time a trajectory inside I takes to 
approach the nonhyperbolic region and return to /. 

n c : is the time when returning without touching the 
nonhyperbolic region becomes as probable as re- 
turning after touching. 

If / is small and away from any nonhyperbolic region we 
expect n* < n a < n c . 

In terms of escapes, the times in Eq. (I25| are inter- 
preted as: 

n*: is a convergence time which is proportional 
to 1/1 A' I, where A' is the negative Lyapunov ex- 
ponent of the saddle (the time to relax to the hy- 
perbolic component of the saddle along its stable 
manifold) . 

n a : is the time needed to approach the nonhyperbolic 
part of the saddle. 

n c : is the time when the nonhyperbolic part of the sad- 
dle becomes as important as the hyperbolic one. 

In this case, the key assumption n* < n a is related to 
the initial density po and leak / being away from the 
sticky regions. Being invariant, the division of the saddle 
is of course independent of the initial density. However, 
for initial conditions touching the KAM island the im- 
portance of the exponential decay can be made arbitrar- 
ily small (see Fig. [8] below) and the power-law exponent 
in (f2"S"]) is modified to a' = a — 1, see Ref. [4g, l5l| (see 
also Ref. Q). 

The above interpretation of p( n) g iven by (|25|) ex- 
presses the view, first suggested in [32j, that the nonhy- 
perbolic component of the saddle is approached through 
the hyperbolic part (32I [HoL HH, [54|, i.e., the fraction of 
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FIG. 6: (Color online) (a) Phase space of the standard 
map J28j with / = [0.25 < x < 0.25 + e, -0.5 < y < -0.5 + e] 
of size e = 0.15, indicated by a checkered box. The KAM is- 
land (nonhyperbolic part of the saddle) is in the center of the 
figure. The support of p r at M(I) is marked as a dark (red) 
region, (b) The unstable manifold of the hyperbolic part of 
the saddle ( black dots) , obtained by using the method of [2(J 
choosing an escape time n s = 80 < n c — 307. The support of 
p r is fully outside the unstable manifold. 

trajectories that reach the nonhyperbolic part is propor- 
tional to those that reach the hyperbolic one. When the 
size of the leak / decreases, the hyperbolic component 
increases while the nonhyperbolic one remains the same, 
meaning that the above picture is even more accurate. In 
this case we can expect that the fraction of trajectories 
that arrive at the nonhyperbolic component is propor- 
tional to the fraction of trajectories that approach the 
hyperbolic component and that therefore the ratio b/a 
[see Eq. ([25]) ] has a weak dependence on 7. Introducing 
this assumption in the definition of the cross-over time n c 
[Eq. ([H] we obtain 

n c ~ I/7 [« 1/(1(1) for small fi(I)], (27) 

meaning that n c is proportional to the reciprocal of the 
escape rate 7. The proportionality constant depends on 
the size of the nonhyperbolic regions. We have verified 
relation (|27[) for different maps, what emphasizes once 
more that for small 7 a well defined exponential decay 
exists for intermediate times n* < n < n c . 

B. Numerical results for the standard map 

We illustrate the previous results in the standard 
map [HIH : 

Vn+i = Vn - 0.52sin(27nz; n ), x n+1 = x n + y n+1 . (28) 

The phase space of this map is shown in Fig. [S] Chaotic 
trajectories stick for an algebraic long time to the hier- 
archical border of the KAM island, shown in the center 
of Fig. [6ja), that constitutes the nonhyperbolic compo- 
nent of the saddle. The exponential decay is governed 
by the hyperbolic component of the saddle, whose unsta- 
ble manifold is shown in Fig. Ob) for a specific leak /. 
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FIG. 7: (Color online) (a) Escape time distributions in the 
standard map with e = 0.15, as shown in Fig. [6] Results 
are shown for po = p r (squares/below) and for po = p M in 
\x\ > 0.25. Lower inset: scaling of n c y with 7 obtained by 
changing e in 0.03 < e < 0.4. Upper inset: log-log plot of the 
main graph over a long period of time, (b) p e , r multiplied by 
exp(7n), where 7 = 0.027403. 



This unstable manifold corresponds to the support of the 
conditionally invariant density p c . Notice its apparent fil- 
amentary character and that it is present inside J, but 
does not enter M(I). 

The distribution of escape and recurrence time for the 
system presented in Fig. [6] is shown in Fig. [Jj As in 
the case of the Henon map (Fig. [3]), short time oscilla- 
tions and an exponential decay are clearly visible. How- 
ever, in this case a slower decay is observed for times 
n > n c w 307. In the upper inset of Fig. [TJa), a log- log 
plot of the main graph shows that the long-time decay 
is roughly a power-law with a w 2.5. We have also ob- 
tained numerically the value of the cross-over time n c 
for leak/recurrence regions I of different sizes by chang- 
ing e defined in Fig. The results are shown in the 
lower inset of Fig.[7[a), and support the theoretical scal- 
ing ([27)) . Remarkably, the value of n c is almost identical 
for recurrence and escape times, the same being valid also 
for other initial distributions po away form the KAM is- 
land. This provides additional support to the picture 
that po quickly converges to the hyperbolic component 
of the saddle, and that the nonhyperbolic component is 
approached in an universal way through the hyperbolic 
component. The transition between the hyperbolic and 
nonhyperbolic components of the saddle is illustrated in 
Fig.il 



VII. FITTING PROCEDURE TO OBTAIN THE 
ESCAPE RATE 

In the previous sections we have seen that the escape 
rate 7 can be formally written as a function of the c- 
measure fi c (I) of the leak region / through relation ([T4"]) . 
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FIG. 8: The saddle characterizing the dynamics at different 
times, ris, given on top of the panels, obtained by the method 
of Ref. |2(J ) . For n s < n c = 307 in (a), the hyperbolic compo- 
nent obtained is clearly disjoint from he region of the KAM 
island. By increasing n 3 as in (b) and (c) the nonhyperbolic 
component appears with an increasing weight. 



However, this relation is of little practical use if one wants 
to determine 7. The reason is that usually only the 
Lebesgue measure is known a priori and the c-measure 
jj, c may change dramatically depending on the position 
of the leak. Its numerical calculation is typically much 
more involved than the determination of 7 itself. In this 
section we search for an efficient method to determine 7 
that goes beyond the rough naive estimate — ln(l— /x(J)), 
given by relations (|13[) and ([6]). 

In general open hyperbolic systems a dense set of 
periodic orbits exists embedded into the nonattracting 
chaotic set. This can be used to effectively obtain 7 

as mm 



-717 



=EIIi A ( r *.»)i' 

r 4 , n 



(29) 



where the product is over the expanding eigenvalues A of 
each periodic point Ti >n belonging to the same periodic 
orbit of length n. The sum runs over all periodic orbits 
of period n that remain inside the system, i.e., outside 
the leak region /. 

We want to take further advantage here of the possi- 
bility of closing the leak of the system. Note that in the 
closed system the decay exponent is 7 = and the sum 
in Eq. ([29| runs over all periodic orbits (inside and out- 
side the leak region). We can therefore give 7 in terms 
of the periodic orbits inside the leak region as 



1 



e ni A ( r ^)i' 



(30) 



From this expression, it is clear that short time periodic 
orbits have a strong influence on 7 since not only the 
primitive orbits, but also all their multiples appear in 
the sum. This dependence is apparent in Fig. 0(b) and 
was previously reported in Refs. [1, [13, The appli- 
cation of ([3T)]) requires, however, again the identification 
of the unstable periodic orbits of high periods. Since the 
deviations of 7 from the estimate 7* = — ln(l — /i(/)) 
has been related to short time periodic orbits we can ex- 
pect to determine 7 from these orbits. A method for 
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this purpose was described in Ref. [27[ in the context of 
recurrences. We adapt it here in the context of escape. 

The distribution of escape times is written as in 
Eq. ©. The oscillatory part for times n < n* will be 
denoted by po(n) and the remaining exponential part 
by Pexp(n) = gexp (— 771). The two parameters g and 
7 of the asymptotic exponential decay can be obtained 
then by imposing the following two conditions: normal- 
ization 



n*-l 



n=l n 

and the value of the mean 



(n) = X np (n) + ^ np exp (n) = 



n=l 



(31) 



(32) 



For recurrences, or p = p r , we can use Kac's lemma (U) 
and thus /ii(I) = p>(I), the natural measure of the leak. 
On the other hand, when initial densities are propor- 
tional to the natural measure in the escape problem, the 
best procedure is to replace (@| by (fl6j) and therefore 
use — Mc(-f) = 1 — e~ 7 on the right hand side of 

Eq. (|3"2"|) . In the last case, we can solve Eqs. (f3"Tj) and 
in terms of 7, which is expressed as the solution of 



(1-Si) 



n*e 7 + 1 - 
eT - 1 



1 



1 - e-T 



(33) 



where Si = S^ =1 1 po(") and 5 2 = J2"=i np (n). Both 
fitting procedures, the one given by Eqs. (|3ip - ([32|) and 
the one using Eq. ([33]) . were employed for the case of 
the Henon map and are marked as dotted lines in Fig. [3] 
Notice that the fitting obtained using Eq. (j3"3"|) alone is 
worse than the one obtained using as additional infor- 
mation the numerical value of (n e ) in Eqs (|3"Tj) - (|3"2"|) . but 
better than the random estimate ©. 

The main advantage of this procedure is that, even if 
theoretically large n* in Eqs. (|31 |) - (f33|) would render bet- 
ter results, in practice the value n* is fairly small as nu- 
merical experience shows (n* (n)) The values of pin) 
for n < n* are thus obtained without much computa- 
tional effort already with a high precision but allow the 
computation of the asymptotic decay 7. Moreover, this 
shows that important information about the long-time 
dynamics is already contained in the short time dynam- 
ics. 

In the nonhyperbolic case one has to take into account 
also the contribution of the power-law tail. In the con- 
tinuous limit, the following contributions: 



^4 = E, C 



7 &( 7 n)- Q = b( 1 n a ) 1 - a /{a-l), 



7716(771) a = 67 



l-a^2-a 



/(«-2), 



(34) 



are to be added as new terms to the left hand side of 
Eqs. (j3T|) and (j32"]) . respectively. Since n a ~ n c ~ l/p(I), 
and a > 2, one sees that these terms become negligible 
for small recurrence/leak regions. 



VIII. CONCLUSIONS 

Experimental and observational measurements often 
occur through holes or leaks that naturally exist or are 
deliberately introduced in an otherwise closed dynami- 
cal system 0, 0, 1 1, 0, 0, 0, H i, G3- The relevant 
observable quantity in such systems with leaks is the 
distribution of escape times from inside the system. In 
strongly chaotic systems, this distributions decays expo- 
nentially. Applying the ergodic theory of transient chaos 
[l| to chaotic systems with leaks an expression for the 
escape rate 7 e [Eq. (fT4|) ] is obtained [5fl in terms of the 
c-measure (associated to the invariant chaotic saddle) of 
the leak I. We have provided a theoretical framework 
to understand the strong dependence of 7 and (n e ) on 
the leak size (Fig. HJ and position (Fig. [5|). In the (un- 
realistic) limit of small leak, the expression of the escape 
rate converges to a previously known relation [Eqs. ©], 
based on the invariant measure of the closed system. Al- 
together our results help to understand previous obser- 
vations @, 0, H3, Hil and set a theoretical framework for 
future applications. 

We have shown that the classical problem of Poincare 
recurrences in closed systems can be described as a prob- 
lem of escape from a system with a leak, once the re- 
currence region is identified with the leak / and the ini- 
tial density is chosen properly. This allows us to treat 
both problems in an unified framework and compare pre- 
vious similar results that were reported independently 
in both fields (compare Refs. [1,13] and |13|), and to 
adapt (in Sec. IVIIj) a previous method to efficiently ob- 
tain 7 e . More surprisingly, we show that the relaxation 
rate of the distribution of Poincare recurrences j r is equal 
to 7 e [Eq. (Tl9|) ] but that the mean recurrence time (n r ) 
given by Kac's lemma Q differs from the typical mean 
escape time [the analog of Kac's lemma for leaked sys- 
tems is (|15j)]. These results are summarized in Table U 
and provide a more detailed account of the results pub- 
lished in Ref. 

In weakly chaotic systems the asymptotic exponential 
decay is replaced by a power-law decay. We have shown, 
however, that if the leak region and the initial distribu- 
tion are away from nonhyperbolic regions, an exponential 
decay is still well defined for intermediate times. We have 
obtained the scaling of the cross-over time n c between the 
two regimes [Eq. (|2"T)) ] which confirms the importance of 
the exponential decay for small leaks. Altogether, these 
results justify an effective splitting of the chaotic saddle 
in a nonhyperbolic component, related to the asymptotic 
power-law decay, and an hyperbolic component, related 
to the intermediate exponential decay for which our pre- 
vious results apply. This splitting is particularly impor- 
tant in the case of Hamiltonian systems (found, e.g., in 
fluid dynamics or optical applications) where the generic 
phase space shows a mixture of regions of regular and 
chaotic motion. 

Finally, it is worth considering briefly two extensions of 
the previous results that have potentially interesting ap- 
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FIG. 9: (Color online) Distribution of the m-th recurrence 
times (m = 1, ...,4, from bottom to top) in the Henon map 
with parameters as in Fig. [3] The solid lines correspond to 
the binomial distribution (|35p using = u c (I). For long 
times, the same decay rate 7 = ln(l — p c (-0) i s observed in 
all curves. The displacement between the distribution and 
the simulation, as in the m = 1 case discussed previously, 
can be related to the short time oscillations that modify the 
pre-factors of the distribution. 

plications. First, consider the case when the leak consists 
of two or more disjoint components (of relevance, e.g., for 
the problem of resetting in hydrodynamical flows @). All 
our results apply to such cases as well, with the remark 
that the escape rate is not the mere sum of the escape 
rates characterizing the components, due to the overlap 
among pre-images, as pointed out by Buljan and Paar [5[ 
and Bunimovich and Dettmann [16| . A second extension 
is to consider the distribution of the second, third, . . . , 
m-th recurrence times. In the notation of Sec. UH the 
m-th recurrence time is defined as T^™^ = rii — rii- r . In 
view of the analogy to escape presented in Sec. [IV] this 
corresponds to the distribution of escape times taking 
place only after entering the leak / for the m-th time. 
The naive (binomial) approximation ([5]) can be extended 
to this case as [27j 

Pr m) ^ = (T (7 V/ )! iM ^ J Hl " M^)) T ~"\ (35) 
(1 — m)!(m — 1)! 



and implies that the asymptotic decay rate 7*( m ) = 
In (1 — [see ((6|)] remains valid for any m. In Fig. [9] 

we show numerical results for the Henon map. For m > 1 
the distributions take, of course, the value zero for n = 1, 
start growing, and then after a maximum is reached, an 
exponential decay sets in. The decay rate has been found 
again to be different from the binomial approximation, 
but apparently also independent of m, i.e., to have the 
value 

7 (m) =7 = Ml-/i c (/)). 

discussed in the context of the first recurrence times [see 
Eq. (fl3|) ]. To understand this result, let us first compare 
the invariant saddles for m = 1 (Si) and m = 2 (S2). 
It is easy to see that Si c £2- Trajectories that belong 
to 5*2 but not to S± must have a point P inside / that 
never returns to I again. Hence, these trajectories nec- 
essarily approach Si for both t — > ±oo. To be outside 
/ for t — + 00, P has to belong to the unstable manifold 
of Si. To be outside / for t — ► —00, the image point 
M[P) taken with respect to the closed system's map has 
to be on the stable manifold of Si. Therefore, points P 
belong to the intersection of the unstable manifold of Si 
and the pre-image (using map M) of the stable manifold 
of Si and also approach S± asymptotically. Accordingly, 
trajectories escaping for long times in the case m = 2 
appear still to be governed by SI having therefore the 
escape rate 7. This argument can be extended straight- 
forwardly to any finite m and can also be important for 
the the case of partial escape in optical systems [l2j . 
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